LIBRARY 



NAVAL POSTGHACUATE SCHOOL 
MONTEREY, CALIFORNIA 93940 



NPS- 59PC7081A 



UNITED STATES 

NAVAL POSTGRADUATE SCHOOL 



/ 



I 




HEAT TRANSFER AND F1DW FRICTION CHARACTERISTICS OF 
SEVERAL PLATE-FIN TYPE COMPACT HEAT EXCHANGER SURFACES 

by 

P. F. Pucci 

// 

S. F. Ball, Jr. 

R. E. Traister 

31 August 1967 



This document has been approved for public release 
and sale; its distribution is unlimited. 



FEDDOCS 
D 208.14/2 
NPS-59PC7081A 



f sf ' 1 

^ jo? -7 2 

' - iT'?/*'’ ' ? /' 



.•r < }.IBP ' 

-Vl. T 
•■}3is4i-o i w 



NAVAL POSTGRADUATE SCHOOL 
Monterey , California 



Captain T. A. Melusky, USN 
Acting Superintendent 



Dr. R. F. Rinehart 
Academic Dean 



ABSTRACT : 

Experimental results for the convective heat transfer and flow 
friction characteristics of plate-fin type heat transfer surfaces 
are presented for five surfaces. Four surfaces had a sinusoidal fin 
geometry. Two of these were fabricated of stainless steel, two of 
nickel. One of the nickel surfaces had perforated nickel fins. The 
fifth surface was fabricated of stainless steel and had a triangular 
fin geometry. The heat transfer data were obtained by the single-blow 
transient testing technique and by the cyclic testing technique. The 
two techniques compliment each other in extending the range of testing 
capability. 
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NOMENCLATURE 



English Letter Symbols 

A Matrix total heat transfer area 

A c Matrix minimum free flow area 

Afr Matrix total frontal area 

A s Matrix solid cross sectional area 

available for thermal conduction 

A^ Matrix solid cross sectional 

conduction area considering 
effect of matrix material 
perforations 

a s Fin thickness 

C c Jet contraction-area ratio for 

circular tube 

Cp Specific heat (gas) at constant 

pressure 

c $ Matrix material specific heat 

Cf Flow stream capacity rate (mCp) 

C s Matrix heat capacity (W s c s ) 

0^ Hydraulic diameter of internal 

passage 

E Friction power expended per unit 

of surface heat transfer area 

G Matrix flow stream mass velocity 

(m/Ac) 

g Proportionality factor in Newton’s 

Second Law 

h Thermal convection surface heat 

transfer coefficient; heat transfer 
power per unit area per degree 
temperature difference 

K c Loss coefficient for sudden flow 

contraction at matrix entrance 



sq ft 
sq ft 
sq ft 
sq ft 

sq ft 
ft 

dimensionless 

Btu/ (lbm deg F) 

Etu/(lbm deg F) 

Btu/ (hr deg F) 

Btu/(deg F) 
ft 

hp/sq ft 
lbm/(hr sq ft) 

32.2 (lbm ft) / (lbf sec 2 ) 
Btu/(hr sq ft deg F) 



dimensionless 
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K 



d 



K 



e 



k 

K 



L 

m 

P 

P 

q 

R 



r 



H 



Momentum velocity-distribution 
coefficient 

Loss coefficient for sidden flow 
expansion at matrix exit 

Unit thermal conductivity 

Matrix material thermal con- 
ductivity 

Thermal conductivity in radial 
direction 

Total matrix flow length 
Mass flow rate 
Pressure 

Matrix porosity (A c /A^) 

Heat transfer rate 

Universal gas constant (53.35 for 
air) 

Hydraulic radius (A C L/A^) 



dimensionless 

dimensionless 

Btu/(hr sq ft deg F/ft) 
Btu/ (hr sq ft deg F/ft) 

Btu/(hr sq ft deg F/ft) 

ft 

lbm/hr 
lbf/sq ft 
dimensionless 
Btu/hr 

(ft lbf)/(lbm deg R) 
ft 



S Relative shape factor 

t Temperature 

u Flow velocity 

V m Matrix volume 

V g Matrix material volume (corrected 

for effects of perforations) 

W Matrix mass 

x Distance from matrix inlet in 

flow direction 



dimensionless 
deg F 
ft/sec 
cu ft 
cu ft 

lbm 

ft 



Subscripts 

atm Local atmosphere 

ave Average 

f Fluid (gas , air) 

i Initial, inlet, individual 
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m 


Mean, matrix 


o 


At orifice 


s 


Solid (matrix material) , static 


STD 


Standard (temperature and pressure) 


X 


Local conditions 


1 


At inlet (upstream of heating elements) 


2 


At matrix inlet (downstream of heating) 


3 


At matrix exit 



Greek Letter Symbols 



P 


Area density: compactness (ratio of 
matrix total heat transfer area to 

core volume-A/V) sq ft/cu ft 


8 


Diameter ratio, do/d dimensionless 


A 


Difference or change (time, distance, 
temperature) 


0 


Time Sec, hr. 


V 


Fluid viscosity (dynamic) lbm/hr ft 


V 


Specific volume ft^/lbm 


0 


Density lbm/cu ft 


0 


Ratio of free flow to frontal area 

(A c /Af ) dimensionless 




Angular frequency of temperature 

oscillation rad /sec 



Dimensionless Groupings 



a 


Amplitude attenuation parameter 


f 


Fanning friction factor (Ratio of wall shear stress to the 
fluid dynamic head). This factor is plotted versus to 

illustrate the flow friction characteristics of a matrix 
surface . 
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j Colburn Factor (N s t N pr ) 5 a generalized heat transfer 

grouping. This factor is plotted versus Np to illustrate 
heat transfer characteristics of a matrix surface. 

M Flow parameter (mc p /Wc s u > 0 ) 

N Nusselt number (hD^/k ), a heat transfer modulus 

Np r Prandtl number (pc p /k ), a fluid properties modulus 

Reynolds number (D^G/p ) , a flow modulus 

N St Stanton number (h/Gcp ) , a heat transfer modulus 

Ntu Number of heat transfer units (hA/mc p ) , a heat transfer 
parameter 

R Ratio of matrix outlet temperature amplitude to matrix inlet 

temperature amplitude during cyclic temperature variation 

X Longitudinal heat conduction parameter for solid matrix 

material (k s A s /mLc p ) . 

5 Distance parameter (x/L) 

Z Time parameter (hA/Ws Cs) used in Single-Blow analysis (17) 

t Time parameter (uO/L ) used in cyclic technique 

u) Frequency parameter (to o L/u ) used in cyclic technique 
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INTRODUCTION 



This document reports the results obtained during FY66 and FY67 
of a continuing program for the evaluation of the heat transfer and 
fluid friction characteristics of compact heat exchanger surfaces at 
the Naval Postgraduate School, Monterey. Results obtained in FY65 
have been previously reported [1], [2]*. 



DESCRIPTION OF SURFACES 

The five surfaces tested, manufactured and provided by Solar, a 
division of International Harvester, San Diego, had basically sin- 
usoidal or triangular fins in a plate-fin geometry. Table I is a 
summary of the matrix geometries. Only one matrix tested, designated 
Solar //3, was brazed. All matrices had a frontal cross section of 
approximately 3 and 1/16 inches square and a length of approximately 
3. inches in the flow direction. The perforated nickel plate fins of 
the Solar #2 matrix is an electro-deposited metallic sheet of integral 
structure manufactured by Perforated Products, Incorporated, and is 
described in Appendix D. 



1. Numbers in brackets refer to numbered items listed in References. 
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EXPERIMENTAL METHODS 



To obtain the heat transfer data, two testing techniques were 
employed: the single-blow, transient testing technique and the cyclic 

testing technique. 

The single-blow, transient testing technique is completely des- 
cribed in reference [2] , however, briefly this method consists of 
heating the test matrix by heated air to a uniform temperature (approx- 
imately 20. F above ambient temperature in these tests) and then sub- 
jecting the matrix to a step change in the air flow temperature to 
ambient temperature. The air temperature downstream of the matrix is 
monitored and recorded versus time. The test rig used in these tests 
is shown in Figures 1-5, and described in detail in Appendix C. A 
series of fine nichrome wire heaters were installed upstream of the 
matrix to heat the air. Turning off the heaters provides the step 
change in the temperature of the air without disturbing the flow. By 
referencing the thermocouples downstream of the matrix to the thermo- 
couples upstream of the heaters, the initial temperature difference can 
be very closely controlled. This temperature difference is continuously 
recorded by a strip chart recorder when a run is made. This recorded 
trace has the distinct advantage that there is no transposition of 
data required which would produce increased uncertainties. Reduction 
of the data follows the method of Locke [6] and Howard [4], whereby the 
Ntu of the surface can be evaluated by determining the maximum slope 
of the fluid temperature difference versus time curve during the cooling 
transient . 



10 



The cyclic testing technique is completely described in reference 
[51. However* briefly this method consists of producing a sinusoidal 
time varying temperature of the steady air flow to the matrix. This 
air inlet temnerature and the time varving temperature of the outlet 
air downstream of the matrix are monitored and simultaneously recorded. 
The test rig used was the same as that used for the above single-blow 
transient tests. A low frequency sinusoidal voltage was superimposed 
on an appropriate d.c. level and fed to the nichrome heaters. Reduction 
of the data follows the method of Bell and Katz [5], whereby the Ntu 
of the surface can be evaluated by determining the ratio of output to 
input amplitudes of the temperature variations. See Appendix A. Again, 
recording these temperatures directly and simultaneously, requiring no 
transposition of data, reduces the data reduction uncertainties. 
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PRESENTATION OF RESULTS 



Heat transfer and flow friction data for each of the matrices 
tested are given in both tabular and graphical form, TABLES I through 
IV and Figures 6 through 15, using the Colbum-j modulus, Fanning friction 
factor, and Reynolds number. 

Figure 16 Is a summary plot comparing the heat transfer and flow 
friction characteristics of the matrices tested. 

Figures 17 and 18 present "figure of merit" type curves. The j/f 
versus Reynolds number thus plotted gives an indication of the matrix 
flow frontal area required for a given pressure drop. 

Figures 19 and 20 oresent curves of heat transfer power versus flow 
friction power on a unit area basis, evaluated for fluid properties at 
standard conditions of dry air at 500. F and one atmosphere. 

CONCLUSIONS 

Of the five matrices tested, the Solar No. 2 matrix with sinusoidal 
fin geometry of perforated nickel, exhibited the best overall performance 
particularly at the higher Reynolds Numbers tested. 

The three stainless steel matrices. Solar Nos. 3, 4, and 5, exhibited 
essentially the same overall characteristics. The triangular fin geometry 
of Solar No. 5 exhibited a slightly higher friction factor with a re- 
sulting lower j/f than the sinusoidal fin geometries. 

The employment of the cyclic testing technique appears to be an 
effective means for extending the range of testing for the facility 
originally designed for the single-blow transient testing technique. 
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Figure 1. Schematic Diagram of Test Apparatus. 
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Figure 2. Test Apparatus with Cyclic Power System. 




Figure 3. Heater and Test Section. 




Figure 4. Test Section with Matrix Removed. 
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Figure 5. Nichrome Wire Heater. 
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SOLAR # 1 



Matrix Material 
Specific Heat , c s 
Thermal Conductivity, k g 
Material Thickness, a s 
Total Heat Transfer Area, A 
Frontal Area, A^ r 
Conduction Area, A s 
Free Flow Area, A c 
Matrix Flow Length, 

Matrix Volume, V m 
Matrix Density, p m 
Compactness , 6 
Porosity, p 
Hydraulic Diameter, 



Nickel 

(Btu/lbm deg F) .106 



(Btu/hr fr deg 
(In) 

(sq ft) 

(sq ft) 

(sq ft) 

(sq ft) 

(ft) 

(cu ft) 

(lbm/cu ft) 

(sq ft/cu ft) 

(ft) 



F) 38.7 
.005 
12.768 
.07022 
.01078 
.05944 
.24917 
.0175 
86.1378 
729.7026 
.8464 
.0046395 



Figure 6. Geometric and Physical Properties 
of Solar #1 Matrix 
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Figure 7. Heat Transfer and Flow Friction Characteristics 
of Solar No. 1 Matrix. 
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SOLAR If 2 



Matrix Material 
Specific Heat , c 
Thermal Conductivity, k g 
Splitter Thickness, a s 
Total Heat Transfer Area 
Frontal Area, Af r 
Conduction Area, A g 
Free Flow Area, A c 
Matrix Flow Length, 1^ 
Matrix Volume , V m 
Matrix Density, p m 
Compactness , 8 
Porosity, p 

Hydraulic Diameter, D^ 
Perforated Fin Material 



(Btu/lbm deg F) 
(Btu/hr ft deg F) 
(in) 

A (sq ft) 

(sq ft) 

(sq ft) 

(sq ft) 

(ft) 

(cu ft) 

(lbm/cu ft) 

(sq ft/cu ft) 

(ft) 

.0047” thick 



Nickel 

.106 

38.7 

.005 

12.093 

.07022 

.01034 

.05987 

.25583 

.01797 

70.3099 

673.143 

.85281 

.0046272 

80/20T 



Figure 8. Geometric and Physical Properties 
of Solar #2 Matrix 
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Figure 9. Heat Transfer and Flow Friction Characteristics 
of Solar No. 2 Matrix. 
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SOLAR it 3 



Matrix Material 


430 Stainless 


Steel 


Specific Heat, c s 


(Btu/lbm deg F) 


.11 


Thermal Conductivity, k 8 


(Btu/hr ft deg F) 


12.8 


Material Thickness, a s 


(In) 


.005 


Total Heat Transfer Area, A 


(sq ft) 


14.3588 


Frontal Area, Af r 


(sq ft) 


.06803 


Conduction Area, Ag 


(sq ft) 


.01325 


Free Flow Area, A c 


(sq ft) 


.05478 


Matrix Flow Length, 1^ 


(ft) 


.24833 


Matrix Volume , V m 


(cu ft) 


.01689 


Matrix Density, p m 


(lbm/cu ft) 


94.2656 


Compactness, 8 


(sq ft/cu ft) 


849.8962 


Porosity, p 




.80517 


Hydraulic Diameter, D^ 


(ft) 


.0037895 



Figure 10. Geometric and Physical Properties 
For Solar it3 Matrix 
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Figure 11. Heat Transfer and Flow Friction Characteristics 
of Solar No. 3 Matrix. 
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SOLAR # 4 



Matrix Material 
Specific Heat, c 
Thermal Conductivity, k 8 
Material Thickness , a 8 
Total Heat Transfer Area, A 
Frontal Area, Af r 
Conduction Area, A 
Free Flow Area, A,. 

Matrix Flow Length, 1^ 
Matrix Volume, V m 
Matrix Density, p m 
Compactness , 8 
Porosity, p 

Hydraulic Diameter, Dj, 



430 stainless steel 



(Btu/lbm deg F) 
(Btu/hr ft deg F) 
(In) 

(sq ft) 

(sq ft) 

(sq ft) 

(sq ft) 

(ft) 

(cu ft) 

(lbm/cu ft) 

(sq ft/cu ft) 

(ft) 



.11 

12.8 

.005 

14.27694 

.06713 

.01229 

.05484 

.24417 

.01638 

88.5918 

P870.969 

.81690 

.00375168 



Figure 12. Geometric and Physical Properti? - 
of Solar #4 Matrix 
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Figure 13. Heat Transfer and Flow Friction Characteristics 
of Solar No. 4 Matrix. 
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SOLAR # 5 



Matrix Material 
Specific Heat, c g 
Thermal Conductivity, k g 
Material Thickness , a s 
Total Heat Transfer Area, A 
Frontal Area, A^ r 
Conduction Area, Ag 
Free Flow Area, A c 
Matrix Flow Length, 1^ 
Matrix Volume , V 
Matrix Density, p m 
Compactness, 8 
Porosity, p 

Hydraulic Diameter, D^ 



430 Stainless steel 



(Btu/lbm deg F) 
(Btu/hr ft deg F) 
(in) 

(sq ft) 

(sq ft) 

(sq ft) 

(sq ft) 

(ft) 

(cu ft) 

(lbm/cu ft) 

(sq ft/cu ft) 

(ft) 



.11 

12.8 

.005 

14.4234 

.06781 

.01148 

.05633 

.26417 

.0179 

81.9345 

805.13313 

.83066 

.0041269 



Figure 14. Geometric and Physical Properties 
of Solar #5 
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Figure 15. Heat Transfer and Flow Friction Characteristics 
of Solar No. 5 Matrix. 
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Figure 16. Comparison of Characteristics. 
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Figure 17. Flow Area Goodness Factors for 
Solar Nos. 1 and 2 Matrices. 




Figure 18. Flow Area Goodness Factors for 
Solar Nos. 3, 4, and 5 Matrices. 
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Figure 19. Heat Transfer Power versus Flow Friction Power 
for Solar Nos. 1 and Matrices. 
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Figure 20. Heat Transfer Power versus Flow Friction Power 
for Solar Nos. 3, 4, and 5 Matrices. 
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Figure 21. NTU versus Amplitude Ratio. 
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Figure 22. Perforated Nickel Geometry. 




Figure 23. Conduction Parameter Analog. 
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TABLE I. Summary Description of Matrices 
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TABLE II. Summary of Friction and Heat Transfer Results 
Solar #1 Matrix 



FRICTION RESULTS HEAT TRANSFER RESULTS 



FRICTION 


REYNOLD'S 


2/3 

N St N Pr 


REYNOLD'S 


NTU 


TESTING 


FACTOR 


NUMBER 


NUMBER 




METHOD 


f 


n r 


.1 


Nr 






.58553 


35.66 


.09482 


35.14 


25.57 


S-B 


.30134 


69.91 


.08829 


68.89 


23.81 


S-B 


.20595 


106.00 


.05464 


104.47 


14.74 


S-B 


.15680 


136.80 


.04155 


134.83 


11.21 


S-B 


.12574 


170.53 


.03475 


168.08 


9.38 


S-B 


.08867 


254.16 


.02215 


250.50 


5.977 


S-B 


.08029 


267.73 


.02033 


263.87 


5.48 


S-B 


.05971 


394.40 


.01401 


388.72 


3.78 


S-B 


.04418 


523.69 


.00959 


516.15 


2.59 


S-B 


.03663 


642.29 


.00775 


633.05 


2.09 


S-B 


.03713 


660.90 


.00812 


660.90 


2.197 


CYC 


.02736 


912.84 


.00525 


899.71 


1.42 


S-B 


.02419 


1060.99 


.00536 


1060.99 


1.449 


CYC 


.02248 


1164.10 


.00394 


1147.36 


1.06 


S-B 


.02064 


1293.89 


.00400 


1275.26 


1.08 


S-B 


.01995 


1346.47 


.00430 


1346.47 


1.162 


CYC 


.01873 


1432.01 


.00380 


1411.41 


1.03 


S-B 


.01773 


1557.13 


.00378 


1534.71 


1.02 


S-B 


.01652 


1676.73 


.00358 


1676.73 


.966 


CYC 


.01681 


1678.22 










.01585 


1854.19 


.00326 


1854.19 


.882 


CYC 


.01486 


2024.35 


.00309 


2024.35 


.778 


CYC 


.01345 


2349.97 


.00279 


2349.97 


.754 


CYC 


.01298 


2504.63 


.00283 


2505.63 


.765 


CYC 


.01265 


2660.83 


.00293 


2660.83 


.792 


CYC 


.01238 


2805.05 


.00299 


2805.05 


.807 


CYC 


.01214 


2945.42 


.00305 


2945.42 


.823 


CYC 


.01202 


3040.39 


.00308 


3040.39 


.833 


CYC 


.01186 


3179.90 


.00310 


3179.90 


.837 


CYC 
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TABLE III. Summary of Friction and Heat Transfer Results 
Solar #2 Matrix 



FRICTION RESULTS HEAT TRANSFER RESULTS 



FRICTION 


REYNOLD'S 


N St N Pr 3 REYNOLD'S NTU 


TESTING 


FACTOR 


NUMBER 


NUMBER 


METHOD 


f 


% 


j % 





.53389 


38.67 


.10087 


38.12 


25.60 


S-B 


.26494 


75.32 


.08623 


74.24 


21.88 


S-B 


.18811 


110.81 


.05455 


109.21 


13.84 


S-B 


.15304 


139.00 


.04330 


136.98 


10.98 


S-B 


.14937 


145.58 


.04233 


143.48 


10.737 


S-B 


.10594 


212.84 


.02602 


209.78 


6.60 


S-B 


.08185 


276.05 


.02066 


272.04 


5.24 


S-B 


.06003 


393.56 


.01317 


393.56 


3.354 


CYC 


.05593 


421.28 


.01217 


415.24 


3.09 


S-B 


.04840 


497.17 


.00926 


490.02 


2.35 


S-B 


.04383 


562.64 


.00966 


562.64 


2.457 


CYC 


.04258 


570.53 


.00679 


562.26 


1.72 


S-B 


.03559 


701.42 


.00801 


701.42 


2.038 


CYC 


.03521 


705.01 


.00598 


694.80 


1.52 


S-B 


.03072 


823.83 


.00513 


812.02 


1.30 


S-B 


.02798 


923.82 


.00641 


923.82 


1.630 


CYC 


.02720 


975.36 


.00428 


961.25 


1.09 


S-B 


.02510 


1064.41 


.00572 


1064.41 


1.453 


CYC 


.03074 


1090.38 


.00531 


1090.38 


1.349 


CYC 


.02505 


1096.74 


.00412 


1096.74 


1.05 


S-B 


.02142 


1368.82 










.02056 


1425.27 


.00535 


1425.27 


1.360 


CYC 


.02078 


1479.07 


.00545 


1479.07 


1.385 


CYC 


.02057 


1538.31 










.02005 


1639.80 


.00602 


1639.80 


1.531 


CYC 


.02021 


1778.11 


.00634 


1778.11 


1.612 


CYC 


.02071 


1887.96 


.00658 


1887.96 


1.672 


CYC 


.02057 


2098.59 


.00671 


2098.59 


1.704 


CYC 


.01848 


2317.65 










.02046 


2318.17 


.00679 


2318.17 


1.726 


CYC 


.02032 


2472.89 


.00688 


2472.89 


1.748 


CYC 


.02019 


2612.20 


.00693 


2612.20 


1.759 


CYC 


.01995 


2821.31 


.00693 


2821.31 


1.760 


CYC 


.01979 


2072.19 


.00690 


2972.19 


1.751 


CYC 


.01951 


3197.52 


.00682 


3197.52 


1.732 


CYC 


.01934 


3355.29 


.00679 


3355.29 


1.724 


CYC 


.01913 


3495.91 


.00674 


3495.91 


1.710 


CYC 


.01890 


3696.46 


.00660 


3696.46 


1.675 


CYC 
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TABLE IV. Summary of Friction and Heat Transfer Results 
Solar #3 Matrix 



FRICTION RESULTS HEAT TRANSFER RESULTS 



FRICTION 


REYNOLD ' S 


M m 2 /3 
N St N Pr 


REYNOLD'S 


NTU 


TESTING 


FACTOR 


NUMBER 


NUMBER 




METHOD 


f 


n r 


j 


H 






.20320 


60.19 


.07076 


59.32 


23.30 


S-B 


. 13761 


91.13 


.04325 


89.37 


14.24 


S-B 


.10327 


121.78 


.03031 


120.03 


9.98 


S-B 


.07288 


172.54 


.01995 


170.06 


6.57 


S-B 


.06223 


221.36 


.01464 


218.18 


4.82 


S-B 


.04422 


313.54 


.01058 


313.54 


3.490 


CYC 


.04132 


342.40 


.00853 


337.48 


2.81 


S-B 


.04015 


357.16 


.00930 


357.16 


3.065 


CYC 


.03507 


401.24 


.00842 


401.24 


2.775 


CYC 


.03125 


456.46 


.00502 


449.90 


1.65 


S-B 


.02705 


531.35 


.00634 


531.35 


2.091 


CYC 


.02574 


567.42 


.00438 


559.26 


1.44 


S-B 


.02427 


601.41 


.00577 


601.41 


1.901 


CYC 


.02200 


672.09 


.00379 


662.44 


1.25 


S-B 


.01938 


809.42 


.00439 


809.42 


1.447 


CYC 


.01737 


902.44 


.00399 


917.52 


1.316 


CYC 


.01421 


1144.49 


.00346 


1144.49 


1.140 


CYC 


.01378 


1194.66 










.01306 


1266.97 


.00331 


1266.97 


1.092 


CYC 


.01229 


1380.58 


.00306 


1380.58 


1.009 


CYC 


.01193 


1436.51 










.01157 


1518.23 


.00285 


1518.23 


.940 


CYC 


.01108 


1675.60 


.00268 


1675.60 


.883 


CYC 


.01123 


1858.27 


.00263 


1858.27 


.866 


CYC 


.00975 


1921.60 










.01121 


2013.75 


.00269 


2013.75 


.886 


CYC 


.01123 


2141.30 


.00285 


2141.30 


.941 


CYC 


.01129 


2220.69 


.00294 


2220.69 


.970 


CYC 


.01134 


2320.36 


.00306 


2320.36 


1.011 


CYC 


.01134 


2487.33 


.00327 


2487.33 


1.080 


CYC 


.01134 


2544.96 


.00330 


2544.96 


1.089 


CYC 


.01136 


2658.22 


.00331 


2658.22 


1.092 


CYC 


.01132 


2807.06 


.00343 


2807.06 


1.130 


CYC 


.01099 


3091.52 


.00342 


3091.52 


1.126 


CYC 



36 



TABLE V. Summary of Friction and Heat Transfer Results 
Solar #4 Matrix 



FRICTION 


RESULTS 


HEAT TRANSFER RESULTS 






FRICTION 


REYNOLD'S 


„ M 2/3 
N St N Pr 


REYNOLD'S 


NTU 


TESTING 


FACTOR 


NUMBER 


NUMBER 




METHOD 


f 


% 


j 


Nr 






.37693 


30.98 


.00611 


30.53 


2.00 


S-B 


.20381 


59.13 


.08412 


58.28 


27.52 


S-B 


.14272 


88.76 


.04310 


87.49 


14.10 


S-B 


.11402 


112.87 


.03061 


111.25 


10.01 


S-B 


.07412 


171.76 


.01845 


169.30 


6.04 


S-B 


.05929 


215.99 


.01499 


212.89 


4.74 


S-B 


.06105 


226.98 


.01370 


223.71 


4.481 


S-B 


.03866 


304.36 


.01139 


306.36 


3.732 


CYC 


.04174 


331.25 


.00847 


326.48 


2.772 


S-B 


.03983 


338.77 


.00806 


333.91 


2.64 


S-B 


.03290 


383.18 


.00939 


383.18 


3.078 


CYC 


.03100 


447.77 


.00467 


441.35 


1.53 


S-B 


.03139 


454.43 


.00468 


447.88 


1.532 


S-B 


.02881 


452.28 


.00826 


452.28 


2.708 


CYC 


.02566 


530.33 


.00735 


530.33 


2.409 


CYC 


.02546 


557.45 


.00429 


549.45 


1.40 


S-B 


.02616 


561.71 


.00700 


561.71 


2.292 


CYC 


.02309 


604.56 


.00671 


604.56 


2.200 


CYC 


.02172 


664.38 


.00346 


654.85 


1.13 


S-B 


.02142 


665.20 


.00625 


665.20 


2.048 


CYC 


.01958 


748.02 


.00577 


748.02 


1.890 


CYC 


.01820 


841.71 


.00535 


841.71 


1.751 


CYC 


.01741 


888.19 










.01756 


891.99 


.00511 


891.99 


1.673 


CYC 


.01441 


1120.75 


.00398 


1120.75 


1.303 


CYC 


.01426 


1162.42 










.01291 


1303.44 


.00358 


1303.44 


1.173 


CYC 


.01239 


1415.34 










.01156 


1526.00. 


.00324 


1526.00 


1.062 


CYC 


.01103 


1668.92 


.00312 


1668.92 


1.020 


CYC 


.01126 


1807.14 


.00310 


1807.14 


1.015 


CYC 


.01116 


1904.64 


.00313 


1904.64 


1.025 


CYC 


.01161 


1908.14 










.01107 


2041.42 


.00321 


2041.42 


1.072 


CYC 


.01106 


2170.84 


.00327 


2170.84 


1.102 


CYC 


.01156 


2557.98 


.00350 


2557.98 


1.144 


CYC 


.01095 


2711.58 


.00351 


2711.58 


1.148 


CYC 


.01088 


2828.75 


.00356 


2828.75 


1.165 


CYC 


.01086 


2887.46 


.00358 


2887.46 


1.172 


CYC 
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TABLE VI. Summary of Friction and Heat Transfer Results 
Solar #5 Matrix 



FRICTION 


RESULTS 


HEAT TRANSFER RESULTS 






FRICTION 


REYNOLD'S 


2/3 

N St N Pr 


REYNOLD'S 


NTU 


TESTING 


FACTOR 


NUMBER 


NUMBER 




METHOD 


f 


% 


j 


Nr 






.22228 


63.05 


.06573 


62.15 


21.15 


S-B 


.14597 


96.54 


.03924 


95.15 


12.62 


S-B 


.11431 


121.54 


.02967 


119.79 


9.54 


S-B 


.08385 


183.46 


.01805 


180.83 


5.81 


S-B 


.06594 


231.11 


.01393 


227.80 


4.48 


S-B 


.04962 


337.08 


.00921 


337.08 


2.972 


CYC 


.04488 


360.72 


.00800 


355.54 


2.57 


S-B 


.03749 


458.12 


.00670 


458.12 


2.160 


CYC 


.03439 


482.55 


.00460 


475.63 


1.48 


S-B 


.03123 


556.72 


.00568 


556.72 


1.831 


CYC 


.02846 


600.21 


.00426 


591.59 


1.37 


S-B 


.02472 


708.74 


.00396 


698.57 


1.28 


S-B 


.02298 


788.74 


.00418 


788.74 


1.347 


CYC 


.01993 


949.26 


.00355 


946.26 


1.145 


CYC 


.01907 


964.34 










.01637 


1203.37 


.00301 


1203.37 


.970 


CYC 


.01548 


1247.65 










.01466 


1405.06 


.00268 


1405.06 


.862 


CYC 


.01341 


1548.85 










.01320 


1635.14 


.00250 


1635.14 


.804 


CYC 


.01304 


1804.55 


.00243 


1804.55 


.783 


CYC 


.01272 


1938.77 


.00247 


1938.77 


.794 


CYC 


.01207 


2040.53 










.01252 


2061.91 


.00249 


2061.91 


.802 


CYC 


.01235 


2214.46 


.00272 


2214.46 


.874 


CYC 


.01233 


2323.94 


.00289 


2323.94 


.932 


CYC 


.01222 


2480.71 


.00319 


2480.71 


1.028 


CYC 


.01218 


2622.09 


.00337 


2622.09 


1.086 


CYC 


.01212 


2736.70 


.00343 


2736.70 


1.104 


CYC 


.01208 


2876.26 


.00353 


2867.26 


1.135 


CYC 


.01204 


2966.69 


.00356 


2966.69 


1.145 


CYC 


.01198 


3028.94 


.00359 


3028.94 


1.155 


CYC 


.01198 


3069.59 


.00361 


3069.59 


1.162 


CYC 
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APPENDIX A 



The cyclic technique employed is that described by Bell and 
Katz (5). The assumptions made in the analysis are: 

(1) All cross sections of the solid normal to some fixed direction 
are uniform* 

(2) The gas velocity is uniform throughout the passage. 

(3) The physical properties of the gas and of the specimen are 
constant. 

(A) The surface heat-transfer coefficient is constant. 

(5) The inlet gas temperature varies sinusoidally. Bell and 
Katz have found solutions to the problem under the following three 
alternative assumptions: 

(a) The thermal conductivity in the solid is infinite. 

(b) The thermal conductivity in the solid is zero in the 

longitudinal (flow) direction and infinite in the transverse direction. 

(c) The thermal conductivity in the solid is zero in the 
longitudinal direction and finite in the transverse direction. 

Of the three solutions, that of case (b) is most nearly satisfied 
in the case of high mass flow rates with its low values of conduction 

parameter. A summary of the theory for this case is presented below: 

Case (b) - Solid having infinite transverse and zero longitudinal 
conductivity . 

Energy balance: 

Fluid: dx + *§q^ = hb(tf - t s ) dx (A-l) 

Solid: P S A S L S II s - dx = hb(t f - t s ) dx (A-2) 

Inlet and Initial Conditions: 
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Inlet : 


t f (0 ,0) 


* t m + A sin co o 0 


(A-3) 


Initial: 


ts(^) = 


tso 


(A- A) 


Initial : 


tf (x,0) 


* to(x) 


(A-5) 



where A is the peak amplitude of the sinusoidally varying tempera* 
ture and t m is some suitable mean temperature. Introducing the 
following dimensionless parameters: 

(1) Time x = ^ 

x 

(2) Position ? “ L 

oj L 

(3) Frequency to = — * — 

% 

(A) Flow: M = m C P 

WsCsoj , , 

hA 

(5) Heat transfer units Ntu ~ ^ 

where : 

u * flow velocity (ft/sec) 

0 = time (sec) 

L = matrix flow length (ft) 

x = distance from matrix inlet in direction of flow (ft) 

u> 4 = angular frequency of temperature oscillation (rad/sec). 

Introducing the above dimensionless parameters into (A-l) through 

(A-5) and rearranging yields: 

Fluid: <5tf Ntu( tf - t s ) = 0 

&K 6 t + 

Solid: 6ts _ M m Ntu(tf — tg) = 0 

6x 

Inlet: tf(0,T) = t so + A sin air 

Initial: t s (0) = t so 
Initial: tf«,0) = t Q U) 



AO 



This system has the periodic solution: 

tf(F,,T) = t m + Ae sin(u)t + QO 

t s (C,T) ■ t m + Aae a ^[M Ntu sin (qjt - tlO - 

cos (urr - fi£)l 



where : 

a = Ntu/ f 1 + (M Ntu) 2 ] 

fl = M(Ntu) 2 /[ 1 + (M Ntu) 2 ] + a) 



The ratio of the inlet temperature amplitude to the exit temperature 
amplitude is 



_ [ Ntu 

R = e 1 1 + (M Ntu) 2 . 



3 



(A- 6) 

This relationshin is plotted in Figure 25 as Ntu versus R for various 
values of M. 

The solutions for the ratio of temoerature amplitudes for cases 
(a) and (c) are presented below: 

Case (a). Infinite conductivity in solid . 

I 

-Ntu. 2 |* 



- e' Ntu ) 2 ]' ; 

R L 1 t H! (1 - e -Ntu) ^ J 



This solution is an approximation to a Taylor series around oj ■ 0 when 
NTU < 2.0. 



Case (c). Solid having finite transverse and zero longitudinal 
conductivity . 

R = exp- 



{ Ntu 

" 1+(M Ntu) 2 



[ , . o 2 M Ntu (M 2 Ntu 2 -l) . H 

1 + 6(1 + m 2 Ntah + - J/ 



where : 

a 2 * 0 ) o as Cs o s /2ky 
as = plate thickness 

ky = thermal conductivity in transverse direction 
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This solution coincides with the value found for R for case (b) 
whenever the series in a can be approximated satisfactorily by its first 
term. For a more detailed explanation on the solutions for the preceding 
three cases the reader is referred to reference (11)* For the treatment 
of the three cases utilizing a liquid as the fluid the paper written by 
Chastain (12) should be consulted. 

On the use of the three solutions. Bell and Katz have shown that 
the ratio of the two amplitude ratios found for cases (a) and (b) are 
coincident when M < 0.1 and NTU 1.0. In the experimental work per- 
formed by Bell and Katz, the ratio of the two amplitude ratios determined 
for cases (a) and (c) was found to be within 1% of unity. Thus, if the 
experimental values of M are selected near 0.1, and the resulting NTU’s 
near unity, no serious error is introduced by assuming zero thermal 
conductivity in the longitudinal (flow) direction, and employing the 
solution for case (b). 
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APPENDIX B 



Data Reduction Relationships 

The relationships presented in this appendix were used to calculate 
the heat transfer and flow friction characteristics of the matrices 
tested. 

Geometry 

The three most significant geometric parameters that provide a 
means of comparison between matrices are so defined that there exists a 
unique relation between them. Determination of any two of the following 
fixes the third: 



(1) Hydraulic diameter 



„ , 4 x free flow area ,Ac „ 

\ " Ar h " . - - - V- — - 4— L 



heat transfer area 



(B— 1) 



(2) Porosity 



P = 



free flow area Ac 



frontal area 
(3) Area Compactness 



(B— 2) 



8 = 



total 

transfer surface area 



A fr L 



total volume 

If (B-2) is divided by (B-3) the previous interdependency results 

r h = f 



(B-3) 



(B-4) 



Maximum Slope 

The maximum slope of the generalized cooling curve is a function 
of NTU and X [A] . 



(I£2^LI±) 

Vt fl ~ 

d (N^) 



- f (Ntu, X) 






where : 



A3 



t = generalized time variable, tj = initial matrix temp. 



T 



hA 

WsCs 



0 



tfi = air temp, upstream of 
matrix, 0 > 0 



Ntu 




tf2 = air temp, downstream of 
matrix 



Therefore : 



1 = ™c p 

Ntu " p 

Furthermore : 

tf2 - H m tfo - tfi 

tfl - ti tfi - t t 



(B-5) 



and : 






(B-6) 



combining (F.-5) and (B-6) 

d / £ r . v - ~- t i | 

V tfi - ti } v. 



sCs 



d hk) 



me. 



_J # d(tf2~tfi) 

tfi — tf dO 



(r.-7) 



considering the following sketch of a generalized cooling curve: 
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d(tf2 - tf) _ j Y I 
d0 " I X | 



X 

chart speed 



d0 



d(tf2 - t fl ) = Y 
(t t - t fl ) = Z 



combining the above measured and recorded quantities with 



WsCs _ matrix heat capacity _ Cs_ -1 
mCp fluid heat capacity Cf 



and equation (B-7) , yields: 

£s_ l_ Y 
C f Z X 



( tf2 - . - H) 

'tfl - HI 



with this value of maximum slope and X 



ksAs 

mc^ 



(3-8) 



enter Table 1 or Figure 3 of Ref, [2] to obtain the corresponding 
value of NTU. 



Cyclic Relationships 

The cyclic tests were made at the higher mass flow rates where 

case (b) of Appendix A is most nearly satisfied. For this case 

D r Ntu 1 

R ■ exp l 1 + (M Ntu) d 



(B-9) 



Solving (B-9) for NTU, we obtain two possible solutions 



NTU = 



1 Wl - (2 H In I/rT? 
2 M2 In 1/R 



(B-10) 



As discussed by Bell and Katz [5] , only the smaller magnitude solution 



has 



physical meaning. 

NTU = - 



Hence 

- >/l - (2M In 1/R) ? 
2 M2 In 1/R 



(B-ll) 



A plot of equation (B-10) is shown in Figure 25. The solid line portion 
of the constant M curves represent equation (B-ll). 
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Mass Flow Rate 



The mass flow rate, m, was calculated in accordance with refer- 
ences (7) and (R) : 



where 



m = 359Kd FaY-/AP ~y (Ibm/hr) 
o o 



(B-12) 



K = flow coefficient, including velocity of approach factor 
C 



= Vl - f L 



with C = orifice coefficient of discharge (20) 



and 3 = ratio of orifice diameter to internal pipe diameter, 
d Q = orifice diameter, inches 

F a = thermal expansion factor for orifice plate 
Y = expansion factor of fluid 



Y = specific weight of fluid flowing = 



RT 



with P = absolute static pressure at orifice ( lb / ft") 
P = 53.3 (ft~lbf/lb -°R) , universal gas constant for air 

i m 

T = [t 0 (°F) + £59.7] °R , abs. temp at orifice 
= pressure drop across the orifice (in. !! o 0) 



making the substitution of K = 



and y ” 



, since the 



Vl - B * 4 ' ?T 

magnitude of the local acceleration of gravity, g, is taken as being 
equal to the magnitude of standard acceleration, gc, into (B-12) yields; 



m = 359 



Vi - p 4 ( 

From reference (7) , Figure 40A 






(FaY) 



(B-13) 



Y = 1 - (0.41 + 0.35 B* 4 ) 



where k = ratio of specific heats for air = 1.4 

= AP p (ln. Ilg abs.) 

P (In. Ilg abs.) 

and from Figure 38 of same reference 
Fa = 1.0 



also : 
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P = (Patm - 13.6 ) .4912 x 144 (lbf/ft^) 



where : 

Patm is the local atmospheric pressure (in. Hg) 

Po is the static pressure upstream of the orifice plate (in. Hg) 
Inserting into (B-13) the required numerical constants to make the equa- 
tion dimensionally consistent and the above defined values or relation- 



ships for P, T, Y, x, k, R and Fa, (B-13) becomes: 
C 



Po 






fn = 589.81 . 

/THF 

Reynold's Number 

By definition: 



d« [1 - (.41 + ^SS 4 ) — 



p APo(Patm - 13.6 



Po 



Patm - 77 t 99.02 

I J • D 



] L t 0 + 459.7 



Nt> = 2b£L 



where G is the mass flow velocity based on the free flow area, Ac 



(B-14) 



(B-15) 



g = t - 9 

Ac 



m 



p Af r 

substituting (B-15) into (B— 16) : 



(B-16) 



there fore : 



N 



m 



R 



uAfr 



Bh 

P 



from (B - 



4) Bl 

P 



1 

8 



Nd =(^-^ ) ( - ) 

R yAfr ’ ( } 

For ease of calculation, substitute (B-3) , 6 = 



Af r L 



into (B-17) 



Np = 



4mL 

Ap 



(B-17) 



(B-18) 
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Fanning Friction Factor 

The following sketch depicts the flow system to be considered: 

la b 2 



YZZZZZZZZZZZZ ZZZ 7 >/T 
c±zz ZZZZZZZZ TZZZZm 
zzzzzz z. zzzzzrzzi jn 
tz^zzzzzzzzzzzzza 
mLZzzzzzzzzzzzzte 



e2zzzzzzzzzzzzzzzi 

ZZZZZZ2 TZZZZZZZZfo, 



la b 2 

For a gas flow heat exchanger, the pressure changes from section 

1 to a and from b to 2 are very small relative to the total pressure, 
therefore the specific volumes v a ^v^, and v^ — V 2 . Since testing is 
performed with air at moderate pressure, the perfect gas law is applicable. 
The relation for the flow stream pressure drop calculation for most 
heat exchanger cores is (10) 



AP 



2gc Pn 



f (K c + 1 - P 2 ) + 2(3. - 1) + f f — 

vi Ac v^ 



entrance 

effect 



flow 
acceleration 

- (1 - P 2 - Ke) 



core friction 






(B-19) 



exit effect 



Solving (B-19) for f In the Isothermal case and recalling from (B-l) 

„ A L _ 1 

that = — also v = — 

f -h-* § - + ff ) h li + ^ i ^ 

where pm = — - - * ■ °2 with the upstream values subscripted by 1 and the 
downstream by 2. Kc is the entrance coefficient and Ke the exit 
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coefficient, both depend on porosity, geometry of the flow cross section, 
and matrix Reynolds' number [9]. Since the ultimate method of calculating 
all values for this report was by digital computer, the analytical 
expressions defined and formulated by Kays in reference [9] were utilized. 

(1) Contraction ratio, C c * 0.611 + .045p + .344 p (B-21) 

(2) Velocity distribution factor, Kd 

(a) circular tubes 
laminar flow 
Kd = 1.333 
turbulent flow 

f = 0.049 (Fanning friction factor, circular tubes) 

Kd = 1.09068 (4f) + 0.0588'JTf + 1 (B-22) 

(b) gap, laminar and turbulent 

Kd (gap) - 1 = 0.6 
Kd (circular) - 1 

(c) square, laminar and turbulent 

Kd (square) - 1 = 1.17 

Kd (circular) - 1 

(d) triangle, laminar and turbulent 

Kd (triangle) - 1 = 1.29 

Kd (circular) - 1 

When utilizing the appropriate Kd for the core under evaluation and for 
the existing type of flow, 

Exit coefficient Ke = 1 - 2Kdp + p^ (B-23) 

Entrance coefficient: Kc = - 2Cc + Cfe(2Kd 1) (b-24) 

Cc 

Considering an order of magnitude approximation, the first term in the 
brackets in (B-20) is the greatest contributor to the friction factor for 
small pressure differentials. 
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Therefore, the approximation that 



pi + p 2 



« PmsPjftPj 



reduces (B-20) to 



the following 

f = [2gcpm || - (Kc + Ke) - || (l+p 2 )] £h 

It is evident from (B-25) that: 

f*2gc pm || p* , 41^ = p/B and G = m/p A fr 
f»2gc pm 

therefore for a given matrix 
fa p 3 / B 

Colburn Factor 

Colbum-j factor is defined as: 

j = Nst NpJ^ = J!_ N pr % 

^ A 

Substituting (B-16) for G and multiplying — yields: 



hA Ac 
mcp A 



Ac » % 

• • Npr^ 3 » 



but NTU 



hA 

mc„ 



therefore : 



j = NTU N pr ^, combining equations (B-l) 



and (B-4) 



Ac = £_ 

A BL 



thus : 



j 



NTU N. 



% 1 



'Pr L 
It is apparent from (B-27) 



£ 

B 



j a 



B 



Heat Transfer Power and Flow Friction Power 

An evaluation of the heat transfer power versus flow friction 
is of interest in that it is a measure of relative performance. 



(B-25) 



(B-26) 



(B-27) 



power 
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Beginning with the definition of the Colbum-1 factor: 



j - N st N pr 

Nst = — *- u - = -_k 1_ 

N r N pr k C p u N r 



therefore : 



solving for h: 



h£h . !_ . N V a 

Cpu N R pr 



D h ' »»i 



(B-28) 



Pr 

If Cp , |j , and N pr are evaluated at standard conditions of dry air 
at 500°F and one atmosphere [10] 
cp = 0.2477 BTU/lb-F 
p = 0.0678 lb /hr- ft 
p = 0.0413 lb/ft 3 
N pr = 0.671 



Equation (B-28) becomes: 



hstd - 0. 02195 (^) (N R j) 



BTU 

hr-ft2 - °F 



(B-29) 



The flow friction power per unit area is derived as follows: 



CO AP , 

f = 2gc pm ^ 


:?) 


rearranging 


AP a.f f 
L 2gc pm 


, since G = 


pmV 


e - AP rh ( v ^ f r 1 


( JL\ , 


f pmVDh \* 1 


E L ^550J- f ptT 


Urn) pT ' 


< U J J 


therefore : 






e -2Tc 


f N3 


(B 


When equation (B-30) 


is evaluated 


at standard conditions (p and 



p m values given above) , the flow friction power per unit area at Standard 
conditions is: 
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e std ■ 1 - 11 “ 10 ' 7 53“ f (?5oo) HP/ftl <B ‘ 31 > 

For comparison purposes the surface geometries were reduced to a 
common hydraulic diameter of 2. x 10 "^ft. 
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APPENDIX C 



Description of Equipment 

The equipment utilized to obtain the data from which the heat 
transfer and flow friction characteristics were determined was basically 
the same for both methods of testing. Only one change was required, and 
that was in the power equipment to produce the necessary temperature 
change corresponding to the method employed. 

This appendix will be divided into two main subsections; 

(1) Equipment common to the cyclic and single blow testing 

systems 

(2) Power generating equipment utilized 

1. Equipment common to the cyclic and single-blow testing systems 
Fluid Source : 

Air, the working fluid, was drawn through the test apparatus by 
a 30HP, multistage, Spencer Turbo-Compressor rated at 550cfm operating 
on a 220 volt a.c. power supply. 

Flow Metering System : 

Flow rates were measured by means of an ASME standard orifice sec- 
tion designed for convenient changing of orifice plates [7]. Thin plate 
concentric orifices with throat diameters from .775 to 2.31 inches were 
employed. The inside pipe diameter, d, was 3.08 inches. Pressure taps 
were located one pipe diameter upstream and one-half pipe diameter down- 
stream from the orifice. 

Fluid Heater System : 

The air heater section consisted of 28 nichrome .0031*' diameter 
wire heaters. Nichrome was selected because of its high resistivity, low 
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thermal conductivity, and low specific heat* 

The nichrome wire was wound on bakelite frames, two heaters to a 
frame. The wires were wound 1/32" apart, 50 and 52 wires respectively 
to each heater, and were formed into a rectangular mesh. The heaters 
on each frame were electrically connected in parallel and then connected 
through a toggle switch to voltage bus bars., 

A schematic wiring diagram for one frame is shown below: 




All heaters selected operate in parallel when power was applied. The 
number of heaters used was dependent on the testing method. For single 
blow testing the number of heaters used was determined by the number re- 
quired to effect a 20°F rise for a given mass flow rate when a 220 volt 
source was applied, e.g., at m = 790. lbm/hr and V » 220 a.c., 22 heaters 
were required. For cyclic testing all 28 heaters were used while the 
input voltage was varied. 

Matrix Holder and Test Section (Figure 4 ) 

The matrix holder and test section were constructed from polyethylene 
plastic. The test section casing contained a removable frame on which a 
set of five thermocouples connected in series was mounted, T 2 . It also 
contained two static pressure taps, one upstream and one downstream from 
the matrix holder. The components making up the casing and holder were 
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machined to close tolerances, thereby minimizing air leakage, and insur- 
ing good alignment of the flow channel through the system. The flow 
channel was 3-3/16" by 3-3/16". Matrices with flow lengths up to 3" 
could be tested. Another removable thermocouple frame containing a set 
of five thermocouples connected in series , , was mounted in the exit 

plane of the matrix holder. The matrices were placed in the holder and 
held in place by insertion of styrofoam plastic on four sides of the 
matrix to insure a snug fit in the holder and proper position relative 
to the gas flow. 

Pressure Measuring System 

Pressure taps were located upstream and downstream from the matrix 
holder and on either side of the orifice. Each tap was connected to an 
appropriate manometer or draft gage via Imperial Company "poly-flow" 
tubing. The following instruments were utilized: 

(1) E. Vernon Hill and Company Type "C" Micromanometer, 0-125" 

(2) Ellison Draft Gage Company, 0-3.0" inclined gage 

(3) Ellison Draft Gage Company, 10" manometer 

(4) Ellison Draft Gage Company, 20" manometer 

(5) Merriman Instrument Company, 50" manometer 

(6) Merriman Instrument Company, 120" manometer 

Any one of the above instruments or combination could be used to measure 
the differential pressures at the orifice or matrix. Cross checking of 
the various instruments assured reliable operation. 

Temperature Measuring System : (Figure 1) 

All temperatures were measured with thermocouples. There were four 
locations in the system where a temperature was monitored; prior to the 
flow straightner, before and after the matrix, and before the orifice. 
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A 30 gage copper-cons tantan thermocouple upstream from the orifice 
plate was referenced to an ice junction. The output was read from a 
Rubicon Company Portable Precision Potentimeter • 

Four sets of 30 gage iron-cons tantan thermocouples were utilized 
at the remaining three positions. Each of the four sets was made up 
with five thermocouples connected in series. All 20 thermocouples were 
of similar construction and were fabricated from the same spools of 
iron and cons tantan wire. Thirty iron-cons tantan couples were made up 
and tested for temperature response (1) with the aid of a Leeds and 
Northrop Company, Type K-3 Universal Potentimeter. From the couples 
tested the 20 thermocouples that had the lease diviation from the average 
output were selected. The maximum deviation from the average output 
of the thermocouples selected was 018 + .2 yv. 

Two of the four sets of iron-cons tantan thermocouples were inserted 
into a 1/4" diameter aluminum tube which served as a radiation shield. 

The tube was mounted in a frame located at the exit of the contracting 
cone. All ten thermocouples in the two sets were individually wrapped 
with teflon tape to prevent shorting. Each set was designated as T^. 
Another set, T^ , of five thermocouples connected in series was securely 
mounted in a frame in the test section located immediately upstream from 
the matrix. The remaining set of thermocouples was mounted in a frame 
located in the exit plane of the matrix holder. The latter set was 
designated T 3 . The output of T^ and T^ were individually referenced to 
the outputs of the two T^'s located upstream at the entrance to the 
system. Therefore, the output of (T^ vs T£) represented the magnitude 
of the temperature rise of the heated air above ambient entering the 
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matrix, and (T^ vs T^) sensed the temperature rise of the exiting gas 
from the matrix above ambient. 

For the single blow tests a record of the output of (T^ vs. T^) 

as a function of time was made for each run. For the cyclic tests the 

outputs of (T| vs T 2 ) and (T^ vs T^) were recorded. The recording 
instrument was a Mosely Company, Model 7100B dual channel strip chart 
recorder. Sixteen fixed full span settings or any magnitude of voltage 
span from 1 millivolt to 100 volts were available for selection. Twelve 

fixed chart speeds in the range of 1 inch per hour to two inches per 

second were provided. 

2 . Power Generating Equipment Utilized . 

For the single blow test runs the power to the heaters was supplied 
by a 240 volt a.c., 60 cycle source regulated by a General Radio Company 
Type W20HM "Variac" auto transformer, 0-240 volt, 8 amp load capacity. 



AUTOTRANSFORMER 




With the arrangement shown above schematically, the amount of volt- 
age supplied to the resistive heaters could be simply regulated. The on- 
off knife switch provided a means by which the power could be quickly 
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secured. 



The cyclic technique required a means of creating a low frequency 
sinusoidal temperature wave. The frequency range was dictated by the 
core material and mass flow rate for each test run. This restriction is 
implied by the requirement of keeping the flow parameter, M < 0.1. At 
the maximum mass flow rate of 2000. lbm/hr used during testing the 
maximum frequency required was 2.1 cycles per second. 

The system finally employed to fulfill the requirements of the cores 
reported consisted of a 1.5 KW, 3 phase, 1725 RPM, inductively wound G.E. 
amplidyne motor generator and a Hewlett Packard, Model 203A low frequency 
function generator. The schematic drawn below indicates the manner in 
which these components were used. 



AMPLIDYNE 




The function generator's voltage, amplitude, and frequency output 
controlled the field excitation of the amplidyne. The output level of 
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the amplidyne was raised above electrical ground by an amount controlled 
by the setting of a high wattage resistance box connected in series with 
a 250 volt coarsely regulated d.c. generator. The magnitude of the d.c. 
voltage fixed the mean temperature, t m , of the flowing gas. 

The frequency response was relatively flat up to approximately 1.2 
cps. Its peak to peak voltage output in this frequency range was 80 
volts. As the system was driven faster by the function generator the 
peak to peak voltage amplitude decreased to about 30 volts at 2.5 cps. 
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APPENDIX D. PERFORATED NICKEL PROPERTIES 



The perforated nickel plate fins of the Solar #2 matrix is an electro- 
deposited metallic sheet of integral structure manufactured by Perforated 
Products, Incorporated, and is designated by the manufacturer as 80/20T. 

The material is .0047 inches thick with 14% open area formed by slotted 
holes nominally .0023 inches by .034 inches, the centers of which are 
spaced .0125 inches in the flow direction and .050 inches transverse. 

For a solid splitter plate or solid fin, the conduction parameter is 
defined as 



where A s is the solid cross-sectional area for conduction which is equal 
to the product of plate width and thickness, and L is the conduction path 
length which is identical to the matrix flow length. The grouping of 
terms r^T represents the thermal conduction resistance, thus, the 
conduction parameter may be expressed as 

1 L 

R th C th k S A S 

Using this form, the conduction parameter for a perforated material 
will be represented by 

x . _! 

k " W k c 

* 

It is useful to define a relative conduction shape factor, S , such 

that 

S* = R th Ak 

( R th) k x 

* 

Thus S is seen to be the ratio of the conduction parameter of a perforated 
material to the conduction parameter of the same material unperforated. 
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By analogy between conduction heat transfer and electric current, 



the ratio of electrical resistances to current flow of two specific 
geometrical shapes, is equal to the ratio of thermal conduction resist- 
ances to heat flow in the same direction. Thus 

L _ Rth 

r k (Rth) k 

Two dimensional models were constructed from electrically con- 
ducting paper (Teledeltos paper) representing repeatable sections of 
the perforated material and its unperforated counterpart. Figure 26 is 
a plan view of the perforated material, the shaded portion representing 
the repeatable section. Figure 27 below it shows the two shapes for 
which the resistances were compared. Measuring the resistances with 
an ohm-meter yielded the ratio of electrical resistances 

— = I*£h = s * = ill 
r k (R t h) k * 

Finally, the conduction parameter for the perforated material is obtained 

X, = S*A where X = £s£§. 

k L C 
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